NONLINEAR SCHRODINGER EQUATION WITH TWO SYMMETRIC POINT 
INTERACTIONS IN ONE DIMENSION 
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Abstract. We consider a time-dcpcndcnt onc-dimcnsional nonlinear Schrodinger equation with 
a symmetric potential double well represented by two Dirac's <5. Among our results we give an 
explicit formula for the integral kernel of the unitary semigroup associated with the linear part of 
the Hamiltonian. Then we establish the corresponding Strichartz-type estimate and we prove local 
existence and uniqueness of the solution to the original nonlinear problem. 
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1. Introduction 

A presence of a certain symmetry in various physical systems is often linked to some symmetric 
double well potential. Models with such kind of potentials are therefore of interest in many fields 
of research. One of the most important phenomena observed in these situations is the spontaneous 
symmetry breaking, which arises in a wide range of physical systems modelled by nonlinear equations 
with double well potential. For instance, in classical physics it has been experimentally observed 
for laser beams in Kerr media and focusing nonlinearity [9] [TO] ■ Spontaneous symmetry breaking 
phenomenon may also arise in the so-called Bose Einstein condensates, where the effective double 
well potential is formed by the combined effect of a parabolic-like trap and a periodic optical lattice 
[6] HU [25] . From a theoretical point of view, nonlinear systems with a double well potential have 
been recently studied in the semiclassical limit [7, 28 or in the limit of large distance between the 
two wells [231 [27]. 

In order to introduce a simplified model of nonlinear Schrodinger equations describing the basic 
features of systems with double-well potentials, we consider in the present paper a one-dimensional 
system where the two wells are represented by two Dirac's S potentials. In other words, we will deal 
with the nonlinear time-dependent Schrodinger equation 

id t ipt = H a ip t + v \ip t \ 2tJ - tp t 

where ipo is the initial data, /i, v are positive numbers and H a , a € R, denotes a one-dimensional 
Schrodinger operator with symmetric delta interactions of strength a and placed at the points x = ±o, 
see section^] below for a precise definition of H a . Although the stationary states and of the symmetry 
breaking bifurcation for equation (jl.ip with cubic nonlinearity have been studied already by Jackson 
and Weinstein [22] . much less is known about the time behaviour of the wave function for fixed 
distance between the two wells. 

The evolution equation (jl.ip with one delta interaction has been considered by several authors 
[2] Q21 [Tni [IZ1 1201 (SI 133] • In 0) among other things, local existence and uniqueness of the 
solution of the associated nonlinear evolution equation is proved. The key point in the analysis of the 
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corresponding nonlinear problem in [2] is a Strichartz type estimate for the time propagator of the 
free Laplacian with one delta interaction restricted to the absolutely continuous part of the spectrum. 
Such a Strichartz type inequality follows, in the model with one delta interaction, easily from the 
(quasi) explicit formula for the time propagator previously given by [3l [4j [18j [31] . 

Here we will deal with the situation when we have two delta interactions of the same strength. As 
in [2] the crucial ingredient is again a suitable Strichartz type estimate. The main technical difficulty 
with respect to the case treated in [2] is the fact that we don't have an explicit expression for the 
associated time propagator at hand. Therefore, as a first step in the analysis of equation we 
derive a formula for the integral kernel of the time propagator of H a in form of a series of certain 
special functions, see Theorem 13. II In the next step we establish the corresponding Strichartz type 
estimate, Theorem 13.51 With the help of this Strichartz estimate and standard tools of the analysis 
of nonlinear Schrodinger equations, [8], we then prove local existence and uniqueness of a solution of 

We would like to point out that our motivation was not only to obtain the decay estimate (|3.3p for 
the evolution operator, but also to give an explicit expression for the integral kernel of the evolution 
operator, which can be used for numerical simulations of equation based on spectral splitting 
methods, see e.g. [29] for numerical analysis of equation with one single Dirac's 6. 

The text is organised as follows. After the preliminary section [5] where we briefly recall some basic 
properties of the operator H a , we announce the main results of our paper, see section [3] Most of the 
proofs are presented in section [4] Some lengthy technical proofs are postponed to the appendices. 
Concerning the notation, we denote by the norm of a function u in L P (R). In the case p = 2 
we will sometimes write ||u|| instead of ||m||2- 

2. Preliminaries 
The operator H a in L 2 (R) is associated to the closed quadratic form 

g[^] = ||^H 2 + a(IV'(a)| 2 + lV'(-«)| 2 ), i>eD(Q) = H\R). (2.1) 
It is well-known, see e.g. [5], that the functions from the domain of H a satisfy the boundary conditions 

ii(±a + 0) = ip(±a - 0) (2.2) 
«V(±a + 0) = V'(±a + 0)-V/(±a-0). (2.3) 
Then H a acts on its domain D(H a ) as 

H a xl) = -ij>", D(H a ) = {ip € H 2 (R\{±a}) : ip satisfies and . 

Notice that, due to (|2.2[) . the function ip e D{H a ) is continuous at x — ±a and therefore D{H a ) 
is a subspace of H 1 (M.). Thus, in the following we denote by ip(±a) the limit (|2.2[) . It is useful to 
recall some basic properties of the spectrum of H a . The essential spectrum of H a is purely absolutely 
continuous and coincides with the positive real axis: 

cress(H a ) = <T &c (H a ) = [0, +oo) . 

From the explicit form of the resolvent operator (H a — see 15, Chap. II. 2] it is straightforward 

to determine the discrete spectrum of H a : if a > 0, then the discrete spectrum of H a is empty. For 
a < the discrete spectrum consists of negative eigenvalues E given given by the implicit equation 

(-2ik + a) 2 = aV 4fca , k = VE, 3fc > . 

In particular, if a < — — then the discrete spectrum of H a consists of only one eigenvalue Ei(a,a) 
defined by 

-Bi (a, a) = [W(-aae aa ) - aaf , (2.4) 
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where W(-) is the Lambert's special function defined by the equation W{x)e w ^ = x, 10J . If a > — ^ 
then the discrete spectrum of H a consists of two eigenvalues: E\(a, a) and 

E 2 (a,a) = -— jj \W(+aae aa ) - aa] 2 . (2.5) 
3. Main results 

We first state a result which gives an explicit formula for the evolution operator associated with 
H a in terms of a series of parabolic cylinder functions. We then make use of this formula to study 
local existence and uniqueness of the nonlinear problem . In order to formulate our results we 
need to introduce some notation. Let P c be the spectral projector of H a on its absolute continuous 
spectrum. Moreover, define 

z n ■= z n (x, y) = Aan+\x + a\ + \y + a\ 

a n := a n (x,y,t) = — , b n := b n {x,y,t) = -j-^= + . 

Theorem 3.1. Let t > and assume that aa ^ — 1. Then the integral kernel of the operator 
exp(—it H a ) P c is given by 

+oo 

U a (t; x, y) = -= e ^~y^ At - — V [p m (t; x, y) + p m (t; -x, -y)] , (3.1) 
V^t 2ir ^ Q 



where 



p m (t;x,y) = (-l) m r^ (t;x,(-l) m y) , 



r„(i; x, y) = s\gn{a)^\a\ 2n+1 (^) e^e^-^D [-2n - 1, -sign(a)2 J (a„ + ib n )} (3.2) 

and D[n, z] denotes the parabolic cylinder function. Moreover, the series on the right hand side of 
(|3.1[) converges for any t > uniformly with respect to x and y. 

Remark 3.2. Here we assume aa ^ —1 for technical reasons. Note that for aa = —1 the second 
eigenvalue E 2 of H a emerges from the continuous spectrum. 

Remark 3.3. The evolution operator e ~ ltH °< then has the kernel 

V a (t; x, y) = <d(-a)e- itEl <p x {y)<pi(x) + 6(-l - aa)e~ ltE ^ (p 2 {y) ^(s) + U a (t; x, y) 

where ifi and ip 2 are the normalised eigenf unctions associated to the eigenvalues E\ and E% and &{•) 
denotes the Heaviside function. 

Remark 3.4. Time evolution in the models involving one or more delta interactions similar to 
ours have been studied, apart from [2], in many other works, see e.g. [Ill 1121 I13j. For example 
in [llj is considered a model of ionization in which to the unperturbed Hamiltonian —-^ri + g5{x), 
with g > is applied a perturbing time dependent potential of the form —grj(t)5(x), where n(t) is 
a periodic function of time. It is then shown that the survival probability of the bound states of the 
unperturbed Hamiltonian tends to zero as t — > oo and the nature of this decay is discussed in detail. 
Similar survival probability is studied in [12j for a model where the perturbing potential is of the form 
rj(t)(S(x + a) — S(x — a)) with r/(t) = rsm(u)t). 

In our model the linear part of the of the Schrddinger equation corresponds to the Hamiltonian H a 
with time-independent potential represented by two symmetric delta interactions. Therefore, instead 
of looking at the time evolution of the (possible) bound states, we focus on the time evolution of the 
" continuous part" of the operator, i.e. on the evolution operator restricted to the absolutely continuous 
part of the spectrum. 
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As a consequence of Theorem 13.11 we obtain the following dispersion estimate 

Corollary 3.5. Assume that aa ^ —1. Then for any there exists a constant C such that for any 
u G L 1 (R) and any t > it holds 



-,-itH a 



\(- 



P c u||oo < Ct-T \\u\\ v (3.3) 



Remark 3.6. The same decay behavior for the model with one 5 interaction was obtained in [2]. 
Dispersion estimates for Schrddinger operators with regular potentials have been well studied in the 
literature, see for example [24j [32] . In [32] it was shown, under certain regularity assumptions on 
the corresponding potential, that the decay rate t~ 1 / 2 of the unitary group e~ ltH as an operator from 
L 1 (R) to L°°(R) is typical for one- dimensional Schrddinger operators. 

Although the assumptions of [24, 32 do not allow to include directly the situation treated in our 
model, it is reasonable to expect that the dispersion estimate (|3.3[) could also be obtained by approx- 
imating the delta interactions by suitable sequences of regular potentials and applying the results of 
[32]. 

Our main result concerning the evolution problem (jl.ll) reads as follows 

Theorem 3.7. Suppose that the assumptions of Theorem \3.1\ are satisfied. Ifipo € H 1 (R), then the 
problem admits a unique local solution 

i) t € C((0,T), H^R)) H ^((O.T), H-^R)) (3.4) 

for some T > 0. Moreover, if ipt is a solution of (jl.ip , then 

HiMI = lhM, «W=*] te(o,r), (3.5) 

where E is the energy functional on H (R) given by 

Sty] = Ua^H 2 + a [\^{a)\ 2 + m-a)\ 2 } + . (3.6) 

4. Proofs of the main results 

We will need the explicit form of the integral kernel of the resolvent of H a . According to [S] 
Chap. II. 2], we have 

([H a -k 2 Y X A{x)= [ K a (x,y;k)<t>(y)dy, <j> e L 2 {R), %k>0, 

where the integral kernel K a is given by 

4 

K a (x, y; k) = K (x, y;k) + J2 K^x, y; k) (4.1) 

with 

Ki(x,y;k) =-(2k((2k + ia) 2 + a 2 e* 4ka )y 1 Li(x,y;k), K {x, y; k) = & e^ x ~^ (4.2) 

and 

Ll(x, y- k) = -a(2k + ia) e **l*+«l e <*li/+«l , L *fa y; k) = L^-x, -y; k) 

L 2 a {x,y-k) = ia 2 e 2lka e i*l«+-l e «lw-»l , L? a (x, y; k) — l? c J K -—x 1 —y\ k). (4.3) 
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4.1. Evolution operator of the linear equation. We split the proof of Theorem 13.11 into several 
lemmas. We start by studying the properties of the function U a (t;x,y) defined by 



U a {t;x,y) = -- I ke- lkH K a (x,y;k)dk. (4.4) 

7T 



Note that 



with 



U a (f,x,y) = U Q (f,x,y) + — [ e' lkH , , /a( ^ ;fc) d k (4.5) 
v ' ,y ' u ^> >^ 2 7t J r {2k + ia) 2 + aV 4fca v ; 



f a (x,y;k) = "S2L J a (x,y;k), U (t;x,y) = J_ 
r-f V47rit 



e i|<E-i/|7« 



Lemma 4.1. For any t > it holds 

! 

2,T 



1 4 

U a {t;x,y) = U {t;x,y) - — ^ Uj(t;x,y), (4.6) 



w/iere 

Proof. By ()4. 5[) we have 

ET.ft *, „) - Ifcfe y) = ± j£ e"^ A 

2^7/ ^ (2k + ia) 2 \ (2k + la ) 2 J ^ 

Since 

f a (x,y;k)=0(k) fc^O, 

for each a;, y, the series on the right hand side of (|4.8|) converges uniformly with respect to k on K. 
We can thus interchange the summation and the integration to get (|4.7j) . □ 

Lemma 4.2. for £ > and aa ^ -1 we have 

oo oo 

Ui(t;x,y) = ^2r n (t;x,y), U 2 (t;x,y) = ^2f n+ i(t;x,y), f n (t;x,y) = -r n (t; x, -y). 

71=0 71=0 

Proof. By Lemma 14. II 

Ux(t;x,y) = Yi 2n+1 a 2n+1 / e^ — __ e **l*-H.| e *fc|»+a| djfe ( 49 ) 

7e ( 2fc + "*) + 



(2fc + ia) 2 



y R v 2Vt 2 ; 

7i—0 \ v / 

oo I — / ■ , \ n 2 

= si 8' n («)y \ a " n+1 (j) e ^ e-( a " +J& ") 2 ^ [-2n - 1, -sign(a)2z(a„ + ib n )] 
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where we apply formula (| A.6|) . see appendix lAl The calculation of U2(t;x,y) follows the same line. 
Let z n := z n {x,y) = z n (x, —y) and define 

a n := a n (x,y,t) = — ^ — and b n := b n [x,y,t) = — ^ + — — . 

As above we obtain 

U 2 (t;x,y) = -p ^+V»+ 2 exp (i (^V) t «H A 2n+2 (^±1, ^ 



n=0 

D 



-2n - 2, -sign(a)2i(a n+ i + ib n+ i) 
where we again applied formula (|A.6[) . □ 
We will also need a uniform estimate on the sequence r n (t; x, y) defined in (J3T2J) . 



r n (f,x,y)\ < Jl |a| 2n+1 f J ) \I n (t;x,y)\, (4.11) 



Lemma 4.3. Let t > be fixed and assume that aa =/= — 1. T/ien i/iere exists a constant C independent 
of t, x and y, such that for all n > 1 we /iawe 

i i 2n n~ 2n 

|r B (t ;a ;,tf)|, |r n+ i(f;z,y)| < C*t* (4.10) 
2 (1 + |a; + o| + \y + a\) M 

Proof. From the definition of a n ,b ni Lemma I A. 11 appendix [AJ and Lemma 14.21 we get 
where 

I n (t;x,y)= ( e- s2 (s + V2a n +tV2b n )- 2n - 1 ds 

Hence 

|i"«(t;x,tf)| < ^\V2b n \- 2n ~ x < C (— 

which implies (|4.10p . The estimate for \r n+ i(t;x,y)\ is completely analogous. □ 

Next we show that U a (t; x, y) defines the kernel of the evolution operator associated to H a . 

Lemma 4.4. Let K a (x,y; k) be the kernel of the resolvent operator of H a and let P c be the spectral 
projection of H a on [0, oo). Then for any test function <p € Cg°(R) we have 

(e- M ^^) (x) = J^U a (t;x,y)0(y)dy. (4.12) 

Proof. Fix <j> G Cg°(R). By 30, Thm.3.1] it holds 

lim e - {lt+e)H ~4> = e- ltH <*4>, (4.13) 

Since the absolute continuous spectrum of H a is the interval [0, oo), the spectral theorem gives 

/•OO 

e -(i t+ e)H a p c = / e -(it+e)X d p^ 
JO 

where P\ is the spectral projector of H a on the interval [0, A]. By the Stone's formula, see, e.g. [261 
ThmVII.13], we have that for any A > 

P([0,A]) = lim / {[H a -(z + iS)}- 1 -[H a -{z-iS)}- 1 }dz. 



5-S-0+ 2ni J 
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Hence 



-(it+6)fe 2 lim /" kK a (x, y; \/k 2 - id) <b(y) dy dk 

-(it+e)k 2 Hm f k Ka U y . yjk 2 + id) (j)(y) dy dk, (4.14) 



1 





e 



1 

71"/ 



where we have substituted k = y/z in the integral containing the kernel K a (x, y; y/ z + id) and 
k = —\[z in the integral containing the kernel K a (x, y; \J z — id), we find out that. Since <j> £ Cq°(R) 
and \kK a (x, y\\Jk 2 + iS)\ is uniformly bounded for S > small enough, we can exchange the limit 
and integration in (|4.14[) . Note also that 

lim kK a ( x, y; y/k 2 ± id) = k KJx,y;±\k\) Vfc^O 

and that kK a (x, y; k) = 0(1) as k — > 0. From the Fubini theorem we then obtain 

f e -(it+e)H ap A (a.) = i. f /" kK a {x,y,k)e' k2{lt+f) dkdy= f U a (t-ie;x,y)(j>(y)dy. (4.15) 

V ' Jl JR JR 

Passing to the limit e — > 0+ in (|4.15|) gives 



iiff„ 



P e ^)(i)= Km U a (t-ie;x,y)<j>(y)dy, (4.16) 



see (|4. 13[) . By (|4. 5[) it can be directly verified that for any x, y £ R it holds 

lim U a (t-ie;x,y) = U a (t;x,y). (4-17) 

£-5-0 + 

Indeed, since Uo(t; x, y) is continuous in t £ C \ {0}, it follows that 

tf«(*;s,l/)-tfa(t-fe; *,») = -- / e- ik2t fa( ?fl ^ dk + o(l) (4.18) 

7r J R {2k + ia) z + a z e l ^ Ka 

as e — > 0+. To show that the first term on the right hand side of (|4.18l) converges to zero, we observe 
that 

lim I* e -*t U*>V>W dk = (4 19) 

J_ R (2k + ia) 2 + a 2 e l4ka v ' ' 

for any R > 0. On the other hand, a direct calculation shows that 
f a {x,y;k) 



a ^ik\x+a\ e ik\y+a\ _|_ g ik\x-a\ £ ik\y-a\ j _|_ QQ.-^ 



(2k + ia) 2 + a 2 e l4ka 

where the term 0(k~ 2 ) depends also on x,y and a. From the asymptotic relation 



00 _ 2 1 / e~ zR \ 

e~ zk - dk — 0\ — — , R -> oo, |arg z| < 7r, 



R 



k \ zR 2 

see e.g. [TJ Chap. 5], we thus deduce that 

zk 2 f a (x,y;k)(l-e- ek2 ) 



R (2k + ia) A + cre MKa 



uniformly with respect to e > 0. The integral corresponding to k £ (— oo, — i?) is treated in the 
same way. This together with (|4.19[) implies (|4.17[) . It follows that for any x £ R the function 
\U a (t — ie\ x, y) \ is uniformly bounded for e nonnegative and small enough and y on compact subsets 
of R. Hence we can interchange the limit with the integral in (|4. 16[) to conclude the proof. □ 



8 HYNEK KOVARIK AND ANDREA SACCHETTI 

Proof of Theorem \3.1\ Since 

U 3 (t;x,y) = U 2 (t;-x,-y), U A {t;x,y) = U 1 (jt;-x,-y), (4.20) 
the statement of the Theorem follows from Lemmata 14. 1[ 14.21 14.31 and 14.41 □ 
4.2. Proof of Corollary 13.51 We need two auxiliary Lemmata. 
Lemma 4.5. Let aa — 1. For any T > there exists a constant Ct such that 

sup \U a {t;x,y)\ < C T t~i, t G [0,T]. 

x,yeV, 

Proof. Clearly 

\U (t;x,y)\< ' 



/47rf 

It thus remains to estimate Uj(t;x,y),j = 1,2,3,4. By (14. 10)) we have 



^2r n+ i(t;x,y) 



< C't*, t G [0,T], 



n=l 

where the constant C depends only on T. On the other hand, for n = it follows directly from the 
definition that 

\ro(t;x,y)\ , \ri(t;x,y)\ < C, 
where C is a real number. □ 

Analogous upper bound on U a (x, y;t) for large t is given in the following Lemma, whose proof we 
postpone to appendix IB"! 

Lemma 4.6. There exists a constant C such that for all t large enough it holds 

sup \Vt U a (x,y;t)\ < C. (4.21) 

x,yGR 

Proof of Corollary \3. 5[ By Lemmata 14.51 and 14.61 we have 

sup \U a (t;x,y)\ < Ct~^ 

x,yGR 

for alH > and some constant C independent of t. Hence 

By density, this inequality extends to all u € L 1 (R). □ 
As a consequence we obtain the following Strichartz-type estimate. 

Corollary 4.7. Let r > 2 and q — -35- Then for any T > there exists a constant C such that 

\\e- UH «P c u\\ LH(0 , T) , Lrm < C\\u\\ 2 VueL 2 (K). (4.22) 

Proof. Let < t < T. In view of Q3.3P and the obvious inequality 

\\e- ltH «P c u\\ 2 < \\u\\ 2 , 

the Riesz-Thorin interpolation theorem implies that 

\\ e -itH a p c M < ct*~* \\u\\ r , - + - = 1 (4.23) 

p r 

for some C andp G [2, 00]. Since H a is self-adjoint, (|4.22l) follows from (|4.23p and [51 Thm. 2.7.1]. □ 
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4.3. Proof of Theorem 13.71 In order to study the well-posedness of the nonlinear equation 
we follow the strategy adopted in [2] . To this end we define the operator T a by 

H a — Ei (a, a) a < 0. 

The key point, apart from the Strichartz estimate, is to prove the following technical result, see also 
[_ Sec.3.7]. 

Proposition 4.8. Let p E [2, oo]. Then the domain D(T a ) of T a is embedded in L P (M.). Moreover, 
there exists a constant C such that for any e < and any u G L P (M) it holds 

IKeTa-lJ-^Hp < C\\u\\ p . (4.24) 

Proof. Consider first the case a > 0. Since i? 1 (R) is continuously embedded in L P (M.) for any 
p e [2,oo], the embedding D{T a ) = D(H a ) ^ L P (R) follows in view of the fact that D(H a ) is a 
subspace of H 1 (R). To prove (|4.24p we note that (eT a — = (e H a — is an integral operator 
with the kernel 

4 

e~~ 1 K a (x,y;iX) = e^ 1 K {x,y]iX) + e" 1 ^ K 3 a (x, y; iX), A = -ie~? > 0. 

3=1 

First we observe that 

e -^, K a W9 )* = ^(e-H,„) W . 

By the Young inequality we thus get 

h' 1 f*Ko(x,y;i\)u(y)dy\\ p < ^ ||e- AH ||iH| p = ||u|| p . 

To control the other terms in K a (x, y; iX), we write 

—a e -A|a:+a| e -A|i/+a| 

^(x,y;*A)+^ ( ^ y . a) = _______ ( 4. 25) 

a 2 e -Mx + a| e -X| B + a | , _ 2A A( | |_| +a |)x 

2A((2A+a) 2 -a 2 e-« c ') V 1 6 e / f ^ U 

a 2 e -A|x+„| e ->| B - a | , A(| +(t |_| d —2 \a\ q< 

see equation (|4.3p . It is easily seen that the function A/((2A+a) 2 — a 2 e _4Aa ) is bounded for A e (0, oo). 
On the other hand, for y < we have < |y — a\ — \y + <z| < 2a. Hence for some constant C > 0, 
independent of y < 0, it holds 

1 _ g-2 An g-A(|i/-o|— |k+o|) i_ g -4Aa 

A >o (2A + a) 2 - a 2 e- 4Aa " A >? (2A + a) 2 - a 2 e - 4Aa " C ' 

The last term in (|4.25p . which corresponds to y > is estimated in the same way. We can thus 
conclude that there exists a constant C QiQ , independent of x,y and A, such that 

\K 1 a (x,y;iX)+K 2 Jx,y:iX)\ < C a , a i e -M*+«l (V A lv+«l + e -A|M-^. (4 26) 

From the Holder and Minkowski inequality we then get 

\\J R (Ki(x,y^X)+Kl(x, y ^X))u(y)dy\\ p <^\\e-^X \\ J^e"^*! + \- a \)\u(y)\ d W ||_ 

< % || e -A|-||| p || e -A|y + a| +e -A|M- a | ||r||u||p 

A 

< ^2 2 +V§ r -* ||tt|| p , (4.27) 
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where | + £ = 1. Since K^(x, y;i\) + K^(x,y;iX) = K^—x, —y; iX) + K^(-x, —y; iX), we arrive at 



n E WlR K i(^y-. iX ) u (y) d y\\p < const hh- ( 4 - 28 ) 

11 3=1 

This completes the proof in the case a > 0. 

For a < we redefine the parameter A in the following way: 

A = -i\Je- x + E 1 (a,a) , A > y/\E^a)\ , 

where Ei(a,a) < is the lowest eigenvalue of H a defined in (12.4[) . Then the integral kernel of 
(sT a — 1) _1 is again equal to e K a (x, y; iX). Moreover, there exists a constant a , such that 



\e- l Kl{x, V ;iX)\ < < ^^^ = e -x\ x+a \ p -x\ y+a \ 



e 



A|e||A- V^oTa;.. 

and similarly for K 3 a {x, y; iX), j — 2, 3, 4. The Holder inequality applied in the same way as in (|4.27l) 
then gives 

oEll/,^,g;<A)"(g)^llp< ™% 



j=l 



Since 



A 2 |e| |A - Vl^iMII > v'PiMI VA > ^^(o, a)|, 
the proof is completed. □ 

Proof of Theorem \3. 7[ The Strichartz estimate (|4.22l) and a straightforward application of the proof 
of [5J Prop. 4.2.1] imply the uniqueness, in H 1 (M.), of a solution to the integral equation 

<pt = e-' itT ° <po-w f e l ^ T « |^| 2 <p s ds. (4.29) 



./o 

Moreover, by Proposition 14.81 and jH Thms 3.7.1] we get the local existence of the solution to (|4.29l) 
and the conservation of energy and charge associated to this equation: 

f[^] + 0(-a)£J L (a J a)|| V » t || 2 = f[ V >o] + e(-a)S 1 (a J a)|| V >o|| 2 J \\<f t \\ 2 = ||<A)|| 2 (4.30) 
for t G (0,T). This immediately implies the local existence and uniqueness of the solution of 

^ = e~ UH ~ Vo - %v [ e l ^ H " \^ s \ 2 ,p s ds (4.31) 



by setting ijj t = e lE i( a > a ) t (p t ii a < and ipt = ft otherwise. Since (|4.31l) is equivalent to (jl.ip . this 
proves the first part of the statement. The conservation of the charge and energy, equation (|3.5[) . is 
a direct consequences of (|4.30p . □ 

Appendix A. 

Lemma A.l. Let n £ N, 7 > 0, 8 ^ 0, and consider the following integral 

A 2n +ih,S)= j e~ ir 2 (r + 7 + ^)- 2ll " 1 dr. (A.l) 



i. IfS> 0: 

A 2n+1 = -i(-2i) n V2^e- l ^ +l5)2/2 D [-2n - 1, (1 - + iS)] (A.2) 
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ii. If 5 < 0: 



A 2n+1 = -i{-2i) n V2^e-^+^ 2 / 2 D [-2n - 1, (1 - i)(7 + **)] 

+z "(ly! e " l(7+4 " )2 ^(- v7(7+ ^) 



(A.3) 



where D(n,z) is the parabolic cylinder function and where H m is the Hermite's polynom of 
deqree m. 



Proof. Let a = —7 — iS, where 5 > 0, and z = r — a, then integral (|A.1|) takes the form 



A _— ia 

A2n+l — e 



~ — iz — 2iza^ — 2n—\ 



dz 



(A.4) 



By means of the change of variable s — e lTr / 4 z and by the Cauchy theorem we then get 

A 2n+1 = e -i«V in7r / 2 / e s2 - 2le " /4aS S - 2 "- 1 ds = e -i« 2 e -W2 / eS 2 -2 4 e" r / 4 a Ss -2«-l ds 



-(e) 



where (£) is the complex path {( e C : £ = (x + £<5)e J7r / 4 , x € K} and where (e) is the complex 
path defined in Figure 19.3 in pQ. Finally, we set q — \/2s and apply formula (19.5.4) by [T] to obtain 



^4-2n+l 



2 (-«)"2" / e 1 2 n-{^e"'*a) qq -2n-l dq 
J{e) 

\-2i) n V2^ie-^ e "' ia ) 2 \-2n - 1, V^e" /4 c 



In the case <5 < then integral (|A.4[) can be written as 



^2n+l = e" 



iz — 2iza„ — 2n—l 



dz 



R-i\S\ 



R+i\S\ 



-iz — 2iza„ — 2n — l 



dz + 2me~ m Res 



e -iz -2i*a z -2n-l jZ = q 



where the integral is computed as the case of 8 > and where the residue is equal to 



Res 



e- ir '{r-a)- 2n -\r = a 



1 d 2 



(2n)\ dr 2n 



__ a (2n)! ds 2 



(2n)! 



' H 2n (Vi< 



□ 



Remark A. 2. By formulas (19.4-6) and (19.13.1) of [T] equations (A.fy) and fA . 3]) can be written 
as 



A 2n +i = -sign(S) i (-2i)"v / 2^e- l(7+4,5)2/2 J D [-2n - 1, s^n(5)(l - i)(y + iS)} 
Remark A.3. If we set 

X = ^— — and T) = 1 fien 



(A.5) 



A 2n+1 = -si 5 n((5) i (-2i) n V2^e -( * +i, » )a £> [-2n - 1, -sign(6)2i(x + if])] . 



(A.6) 
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Appendix B. Proof of lemma I4~6l 
Since \Uo(t; x, y)\ = l/2\/7rt, it suffices to consider Uj(t;x,y), j = 1,...,4. Let 

4 

f a (x, y;k) = L i( x > V'i k ) = 2 ^kq(x, y; k) + ia 2 p(x, y; k), 
j'=i 



with 



q{x, y; k) = _ e ik\x+a\jk\y+a\ + jk\x-a\ & ik\y-a\ ) ^ y . fc) = ^ y . fc) + h (_ x> _ y . fc) 
h(x, y; k) = e <fc l x+a l 



e i2fea e ifc|y-a| £ ik\y+a\ 



and define 



r a {k) = [l + a 



i [In 



- 1 



4«a + 4fc. 



With this notation the statement of the iemma will follow, if we prove 
Lemma B.l. Let 



V^x,y) 



e -ik 2 t ia 2 p(x,y,k) ^ 
kr a (k) 



Then there exists a constant C such that for t large enough 



sup 



\V*(x,y) + W t a (x,y)\ < Ct^ 2 . 



(B.l) 



Proof. By means of the change of variable 



we get 



C = kVt- 



\V*(x,y)\ < 



2V~t' 



C 



z = \x + a\ + \y + a\ 



Vt 



r a ((/Vt + z/2t 



We can estimate the above integral with the help of the Cauchy theorem for analytic functions 
replacing the domain of integration with the path 



7=<meC : w — ( — S 



it 



8 > 0, C € 



Since the function r a (k) has at most two zeros, both of them having imaginary part strictly positive, 
it follows that for t large enough |r Q (w)| is bounded from below by some positive constant in the 
strip \^w\ < S . A direct calculation shows that for t large enough and some Cs < oo we have 



/ 



1 



r a (CM + z/2t) 
which yields the desired estimate for V*: 



< c 



_ 25 C 

e d ( < C s , 



sup 



K(x,y)\ < Cr 1 ' 2 . 



(B.2) 



To prove the analogue of (|L3.2|) for , we first observe that 

h(x, y; k) = if y > a. 
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We therefore define 

hi(x,y;k) = h(x,y;k) = e ^+a\ e ik\ v - a \ 2 i S m(2ka) if y < -a 



h 2 (x,y;k) = h(x,y;k) = e ik \ x+a \ 
and consider separately the integrals 

Ij(x,y,t) -- 

For the first integral we have 



e i2ka e ik\y-a\ _ & ik\y+a\ 



= ke lk \ x+a \S(y-k) ii\y\< 



e-^ h j {x ' y : k) dk, j = 1,2. 



h(x,y,t) 



kr a (k) 



kr a {k) 



dk 



_ lfe 2 t e 4fcz 2isin(2fca) 



kr a (k) 



dk 



Aia 



iz 2 /it / tQ 



y2 sin 



in [2a(C/Vt + z/2t)] 



sft Jr^ 2a((/Vi + z/2t) r a (C/Vt + z/2/t)' 

where z = \x + a| + \y — a\. Note that the for any 8 > there exists a constant Cs such that 

sin [2a(C/Vt + z/2t)] 



< C s |9fCl < & 



2a{C/Vt + z/2t) 

We can thus mimic the argument of the proof of Lemma IB. II to arrive at the estimate 

\h(x,y,t)\<C/Vi, 

which holds for all t large enough uniformly in x and y. The same arguments apply to the integral 
h{x,y,t) - 



e _ ikH h 2 (x,y-,k) dk 



kr a {k) 

J\x+a\ 2 /4t f i( 



e -ikH e ik\ x+a \S(^_ dk 
r a (k) 



S{y;(/Vt+\x + a\/2t) 
r a ((/Vi+\x + a\/2t) 



d( 



since 



\S(y;w)\<C, Vye [-a, a], \$Sw\ < 8. 
Hence we can conclude that for t large enough 

sup K(x,y)| < CT 1 ' 2 . 



(B.3) 

□ 
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